Isoscaling behavior in the Fission Dynamics 
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The fission processes of "^sn + "2sn and "«Sn + Sn are simulated with the combination of 
the Langevin equation and the statistical decay model. The masses of two fission fragments are given 
by assuming the process of symmetric fission or asymmetric fission by the Monte Carlo sampling 
with the Gaussian probability distribution. From the analysis to the isotopic/isotonic ratios of the 
fission fragments from both reactions, the isoscaling behavior has been observed and investigated 
in details. Isoscaling parameters a and /3 are extracted as a function of the charge number and 
neutron number, respectively, in different width of the sampling Gaussian probability distribution. 
It seems that a is sensitive to the width of fission probability distribution of the mass asymmetrical 
parameter but /3 is not. Both a and /3 drop with the increasing of beam energy and the reduced 
friction parameter. 

PACS numbers: 24.75.+i, 25.85.Ge, 21.10.Tg 



I. INTRODUCTION 

The availability of exotic nuclear beam with extreme 
neutron-to-proton ratio provides an opportunity to ex- 
plore the colhsion dynamics of isospin-asymmetric nu- 
clear systems Q. To facilitate this kind of study, the 
suitable selection of the sensitive experimental observ- 
ables related to the isospin degree of freedom is one of key 
points. One of such observables is the isotopic/isobaric 
ratio 0, 13, which has been used to probe the isospin 
equilibration at medium energies before. Recently, this 
kind of ratios have been systematically revisited for two 
different reactions with the same charge number and the 
similar temperature and a so-called isoscaling law has 
been observed experimentally 0, 1511, ■ Isoscaling means 
that the ratio of isotope yields from two different reac- 
tions, 1 and 2, i?2i(iV, Z) = Y2{N, Z)/Yi{N, Z), is found 
to exhibit an exponential relationship as a function of the 
neutron number N and proton number Z Q 



R2i{N, Z) 



Y2{N, Z) 
Yx{N,Z) 



CexpiaN + I3Z), (1) 



where C, a and P are three parameters. In grand- 
canonical limit, a — A/i„/r and /? — A/i^/T where 
A/i„ and A/i^ are the differences between the neutron 
and proton chemical potentials for two reactions, re- 
spectively. This behavior is attributed to the differ- 
ence of two reaction systems with different isospin asym- 
metry. It is potential to probe the isospin dependent 
nuclear equation of state by the studies of isoscaling 
0. So far, the isoscaling behavior has been experimen- 
tally explored in various reaction mechanisms, ranging 
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from the evaporation Q , fission and deep inelas- 
tic reaction at low energies to the projectile fragmen- 
tation 1 1( 1 1 111 and multi-fragmentation at intermediate 
energy |4L ll2 . lli^ . While, the isoscaling phenomenon has 
been extensively examined in different theoretical frame- 
works, ranging from dynamical transport models, such 
as Blotzmann-Uehling-Uhlenbeck model Quantum 
Molecular Dynamics model 1 141 and Anti-symmetrical 
Molecular Dynamics model [13, to the statistical mod- 
els, such as the expansion emission source model, the 
statistical multi-fragmentation model and the lattice gas 
model ^, Ji,^, JJZ.,.18J. In this work, we will focus on the 
detailed simulation studies on the isoscaling behavior of 
the fission fragments. A brief report has been published 
recently 

In this work, we present an analysis for the fragments 
from the fission which was simulated by the Langevin 
equation. The isotopic or isotonic ratios of the different 
fragment yields from "'^Sn -I- "^Sn and "^Sn -I- "^Sn 
system are presented and the features of isoscaling be- 
havior in fission dynamics are investigated. 

The paper is organized by the following structure. In 
Sec. II, a brief description of the Langevin model is given 
and the partition of masses of two fission fragments is 
assumed; In Sec. Ill, the detailed results for the fission- 
fragment isotopic and isotonic distribution are presented 
and the isoscaling behavior is explored; Finally we sum- 
marize the present work. 



II. BRIEF DESCRIPTION OF THE LANGEVIN 
MODEL 

The process of fission can be described in terms of col- 
lective motion using the transport theory 20, 21, 22, 2 J, 

. The dynamics of the collective degrees of freedom is 
typically described using the Langevin or Fokker-Planck 



2 



equation. In this work, we deal with a Combined Dynam- 
ical and Statistical Model (CDSM) which is a combina- 
tion of a dynamical Langevin equation and a statistical 
model to describe the fission process of heavy ion reac- 
tion 21]. This model is an overdamped Langevin equa- 
tion coupled with a Monte Carlo procedure allowing for 
the discrete emission of light particles. It switches over to 
statistical model when the dynamical description reaches 
a quasi-stationary regime. We first specify the entrance 
channel through which a compound nucleus is formed, 



i.e. the target and projectile is complete fusion. The fu- 
sion process of simulating the fission in each trajectory 
with angular momentum L — til \s described by 



a{l) 



27r 



2/ + 1 



fc2 1 + exp[{l - QI5i] 



(2) 



where the parameters Ic and 51 is according to an approx- 
imating scaling of Ref . |2^ . Namely, 



J 



A. 



1/3n 



(0.33 + 0.205 X V-E^c.m. - K) 



(3) 



r 



when < i?c.m. - V'c < 120 MeV; and when Sc.m. - K > 
120 MeV the term in the last bracket is put equal to 2.5. 
In the above equation, At and Ap represents the mass of 
target and projectile, respectively, and Aqn is the mass 
of compound nucleus. For the barrier the following 
ansatz is used: 



7C3 X 



A 



1/3 



A. 



1/3 



1.6 



(4) 



with C3 = 0.7053 MeV. The diffuseness 51 is found to 
scale as 



51 



{ApAr f'^ X 10-5 X [1.5 + 0.02 x {E^,^ 



^ {ApATf2 



X IQ-^ X [1.5 - 0.04 X {Ec 



Vc 



10)] for 
10)] for 



< K 



10 
10 



Trajectory with the particular angular momentum L 
is started at the ground state position q^s of the entropy 
S{qgs, -Etoti Z,L), q is half of the distance between the 
centers of masses of the future fission fragments. In this 
work the total initial excitation energy E^^^ is given by 
= EbcamAT/{AT + Ap) + Q whcrc Q is the fusion Q- 
value calculated by Q = Mr+Mp-M^^. Mt and Mp is 
the mass of projectile and target come from experimental 
data, respectively. If it is unavailable, it is calculated by 
macroscopic- microscopic model is the mass 

of the compound nucleus which is calculated from the 
liquid-drop model. 

The dynamical part of CDSM model is described by 
the Langevin equation which is driven by the free energy 
F. F is related to the level density parameter a{q) |2y] 



Fiq,T)^V{q)~a{q)T^ 



(5) 



in the Fermi gas model, where V{q) is the fission potential 
and T is the nuclear temperature. 

The overdamped Langevin equation reads 



dq 
'di 



1 



,dFiq,T)T. 



VDiq)m, (6) 



where q is the dimensionless fission coordinate defined 
as above. Po{q) is the reduced friction parameter which 
is the only parameter of this model. The fluctuation 
strength coefficient D{q) can be expressed according to 
the fluctuation-dissipation theorem: 



D{q) 



T 



(7) 



where Af is the inertia parameter which drops out of the 
overdamped equation. T(t) is a time-dependent stochas- 
tic variable with Gaussian distribution. Its average and 
correlation function is written as 



< r(t) >= 0, 

< r(t)r(<') >= 2S,{t-t'). 



(8) 



The potential energy V{A, Z, L. q) is obtained from the 
finite-range liquid drop model [27| 



V{A,Z,L,q) = a2[l-k{ 



Z 



N - 



A 



]A^/^[Bs{q)-l] 



+C3j^[B,{q) - 1] + CrL^A~^/^Br{q), (9) 



3 



where Bs{q), Bc{q) and Br{q) means surface, Coulomb 
and rotational energy terms, respectively, which depends 
on the deformation coordinate q. 02, C3, A: and Cr are 
parameters not related to q. In our calculation we take 
them according to Ref. pClj|. 



a2 = 17.9439MeV, cg = 0.7053MeV, 
k = 1.7826, Cr = 34.50MeV. 

We use c and h to describe the shape of nucleus, 
p\z) = (1 - - ^)c2 + S3h(c, h)z^). (10) 



where 



Co 



= cR, R = I.IQA^'^. 



(11) 



The nuclear shape function Bsh (c, h) and the collective 
fission coordinate (/(c, h) of mass number A is expressed 
as 



Bsh(c, h)^2h 



c- 1 



g(c,/i) = ^c(l + ^Bsh(c,/j)c3). (12) 



The fission process of the Langevin equation is propa- 
gated using an interpretation of Smoluchowski 32] (A — 1 
in the following equation) which is consistent with the ki- 
netic form which reads 



qn+i ^ qn + [ 



-\nT + A[— ( N ..f )]«'^ 



Mq)M dq 



'dq'Mq)M' 



r 



T{q) . 
'Po{q)M^ 



(13) 



Here r is the time step of the Langevin equation, w„ 
is a Gaussian distributed random number with variance 
2. S{q) = 2^a(g, A) [Etot - V{q, A, Z, I)] is the entropy. 
The parameter A allows us to distinguish between the 
different possibilities to discretize the Langevin equation. 
It is called interpretation in the literature. In the anal- 
ysis of the experiments on fission of hot nuclei discussed 
in the review [23, and in the papers quoted there, 
the Ito-interpretation (A = 0) [s^l has been used exclu- 
sively. Also there are other interpretations, namely that 
of Stratonovich jsj (A = 1/2), or an interpretation which 
is consistent with the kinetic form of the Smoluchowski 
equation of [l^ (A = 1). In this work, we take A = 1. 

In our calculation we adopt one-body dissipation 
(OBD) friction form factor /3obd as Poiq) which is 
calculated with one-body dissipation with a reduction of 
wall term except the special case which we claim. Here 
we use an analytical fit formula which was developed in 
Ref. H, i.e. 



. , . . 15/q"-43 + 1 

/3obd('?)- {32_32.2iq 



10.5q^-^ + q^ iiq> 0.38 
if g < 0.38 



In the dynamical part of the model the emission of light 
particles (n, p, d, a) and giant dipole 7 are calculated at 
each Langevin time step r, the widths for particle and 
giant dipole 7 decay are given by the parametrization of 
Blann js^ and Lynn [s^, respectively. 



III. RESULTS AND DISCUSSIONS 

A. Isotopic/Isotonic distributions of the fission 
fragments 

In Fig. n (a) and (b), we demonstrate that the ratio 
of the pre-scission neutron number R^ and of the pre- 
scission proton number Rp between ^^^Sn + ^^^Sn and 
^^^Sn -|- ^""^^Sn, respectively, as a function of beam energy 
{Ehcam/A). First, the values oi Rn is larger than 1 while 
those oi Rp is less than 1, indicating that the neutron 
is easier to be emitted for neutron-rich system while the 
proton is in contrary trend. Of course, this is a natural 
result from the chemical composition of reaction system 
37, 38J. Second, with the increasing beam energy, Rn 
shows a decreasing trend while Rp in reverse way, which 
can be interpreted the isospin effect weakens as the beam 
energy rises up. 

Within the framework of the Langevin simulation we 
chose 200,000 fission events which happen on dynamic 
channel (we give up the events which happen in statis- 
tical part of CDSM model) and chose a Gaussian dis- 
tributed random number as the mass asymmetry param- 
eter ao = aI+aI ^^^^ the system reaches to the scission 
point. When ao = it means the symmetrical fission. 
It is taken from a Gaussian distributed random number 
from -1 to 1 with the mean value of zero. Ai and A2 
is the mass of the two fission fragments, respectively. In 
this work we assume the fission fragments have the same 
N/Z as the one of the initial system and then Zi or Z2 
of fission fragments can be deduced from Ai or ^2- This 
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FIG. 1: The ratios of the pre-scission neutron number (a) and of 
the pre-scission proton number between ^^^Sn + '^^^Sn and ^^■^Sn 
+ -'^-'^Sn as a function of beam energy (i^bcam/^)- 
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FIG. 2: (Color online) The fission-fragment mass distributions 
produced by the Langevin simulation for the reactions of ^^^Sn -|- 
^i^Sn (open symbols) and "^Sn + "^Sn (filled symbols) at 8.4 
MeV/nucleon with the different sampling width ((Tao) of Gaussian 
probability distribution. 



assumption is similar to the case of deep inelastic heavy 
ion collisions at low energies, where the isospin degree 
of freedom has been found to reach equilibrium first "39j . 
Fig.[21shows the mass distribution of the fission fragments 
from ^^^Sn + ^^^Sn and ^^^Sn -|- ^^^Sn reaction systems 
assuming the different width of the sampling Gaussian 
probability for the mass asymmetrical parameter of fis- 
sion fragments (ctqo). Naturally, the bigger the Cqq, the 
wider the fragment mass distributions. 

Samples for the isotopic and isotonic distributions in 
some given Z and N are shown in Fig. |31 The square 
of the full widths of these distributions shows a system- 
atic increase with the Z or as shown in Fig. 0] and the 
absolute value of the differences of the centroid of the 
isotonic/isotopic distributions shows an increasing trend 
too (see Fig. Apparently, the widths are not sen- 
sitive to the width of the Gaussian probability for the 
mass asymmetrical parameter of fission fragments, but 
the differences of the centroid of the isotonic/isotopic dis- 
tributions shows the dependence on it. 

From a practical point of view, the isoscaling behavior 
occurs when two mass distributions for a given Z from 
two processes with different isospin are Gaussian distri- 
butions with the same width but different mean mass. In 
this case, the isotopic distribution in a given Z, namely 
Y{N)\z, and isotonic distribution in a given N, namely 
Y{Z)\j^, can be described by single Gaussian distribu- 




FIG. 3: (Color online) The isotopic (a) and isotonic (b) distribu- 
tions of fission-fragments in some given Z or N (see texts in the 
inserts) from the collisions of ^^®Sn + ^^^Sn (filled symbols) and of 
"2Sn -I- ll2Sn (open symbols) at 8.4 MeV/nucleon. o-qq = 0.06. 
Notice that the scale of X-axis is different. 
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tion, respectively, i.e.: 



Y{N)\z^exp[ 



Y{Z)\n ^ exp[ 



(N - Nz? 
24 

2a^ 



where Nz and iVjv are the centroid of isotopic and iso- 
tonic distributions, a\ and cr^ describe the variance of 
distributions for each element of charge Z and neutron 
number iV, respectively. This leads to an exponential 
behavior of the ratio R21 if the quadratic term in Nz or 
iVjv is neglected, it reads 



ln{R2i{N)\z) 



ln{R2i{Z)\N) 



[{Nz)2 - {Nz)i]N 



[{Nn)2 - {Nn)i\Z 



'N 



(15) 



or tr^ to 



Note that Ea. H15() requires the values for a\ 
be approximately the same for both reactions, which is 
a necessary condition for isoscaling. Indeed, we observed 
this case in our simulations for both Sn + Sn collisions. In 
the Langevin equation, (t| or cr^ essentially depends on 
the physical conditions reached, such as the temperature, 
the density and the friction parameter etc. Considering 
that i?2i(-/V)|z ~ exp{aN) or R2i{Z)\n ^ exp{pz) for a 
given Z ov N, we can get 
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FIG. 4: (Color online) The square of full width of the isotopic (a) 
and isotonic distribution (b) as a function of Z and A^, respectively, 
in different CTcq . The width is the average value of the two systems 
(they are almost the same for both systems). 



a ' 



{Nz)2 - {Nz)i 



{Nn)2 - {Nn)i 



'N 



(16) 



Assuming other ingredients can be neglected, tT| or cr^ 
could be considered to be proportional to temperature 
T of the fission-fragments according to the fluctuation- 
dissipation theorem in this circumstance, 



/3' 



{Nz)2 ~ {Nz)i 
T 

{Nn)2 - {Nn)i 



(17) 



where [{Nz)2 — {Nz)i\ or [{Nm)2 — {Nn)i\ can be under- 
stood as a term of the average difference of the neutron 
or proton chemical potential between two reactions. 

As we already showed that both cr| and ct^ rise with 
Z and TV, respectively, and are almost independent of 
(Jao from Fig. Q^a) and (b) in our calculation. On the 
other hand, we recognize that the similar behavior of the 
Z-dependence of az has been experimentally observed 
in the spallation-fission data of ^"^Pb (1 GeV/nucleon) 
+ d or p etc in Gesellschaft fiir Schwerionenforschung 
(GSI) 40, 41]. According to the model which is based on 
the modern version of Abrasion-ablation model involving 
the fission nuclei by Benlliure et al. j43 |. the square of 



the width of symmetric fission fragment from the macro- 
scopic potential can be expressed by 



4 



1 



Trs 



2 0lCn 



2Cn 



(18) 



where E^^ is the excitation energy above the the fission 
barrier, a is level energy parameter, Tfig is the temper- 
ature of fissioning nuclei and Cmac is the curvature of 
macroscopic potential energy Knac as a function of charge 
asymmetry. In this way, the width of symmetric fis- 
sion fragment distribution increases with temperature. 
This is also the case in our present model calculation. In 
the other word, the temperature of the fission-fragments 
which mostly originates from the symmetric fission seems 
to increases with the charge number of fragments. Re- 
cently, a systematic study on the experimental data also 
displays that the variance of the fragment mass distri- 
bution increases with the temperature of the compound 
nucleus and the fission- fragments 43]. 

To verify the relationship of the temperature and 
charge number of fragments as stated above, we extract 
the temperature of the fissioning nuclei in the scission 
point when the system happens on dynamic channel. 
Fig. Ufa) and (b) demonstrate that the mean temper- 
ature of two systems as a function of Z and N for the fis- 
sioning nuclei, respectively. Obviously, the temperature 
almost increases linearly with the charge number (Zgs) 
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FIG. 5: (Color online) The extracted mean temperature of the 
fissioning nuclei of two reaction systems as a function of their charge 
number (top panel) or the neutron number (bottom panel). 



FIG. 6: The absolute value of {Nz)2 - {Nz)i (a) and \{Nn)2 - 
{^n)i\ (t>) as a function of Z and of the fission- fragments, re- 
spectively, between ^^^Sn -|- ^^''Sn and ^^^Sn + ^^^Sn with different 



or neutron number (-/Vfis) of the fissioning nuclei. Since 
we assume the fission-fragments have the same N/Z as 
the one of the fissioning nuclei, hence the temperature 
of the fission-fragments shall increase with their charge 
number. 

In Eq.(16) [iNz)2 - iNz)i] or [(7V^)2 - (iVjv)i] can 
be understood as a term of the average difference of 
the neutron or proton chemical potential between two 
reactions. Fig. E^a) and (b) shows [{Nz)2 — (^z)i] 
and the absolute value of [{N]y)2 — {Nn)i] as a func- 
tion of Z or N, respectively, in different ■ Ap- 
parent increasing behavior with Z or N has been ob- 
served. In order to understand the increasing behavior 
of [{Nz)2 - iNz)i] or \[{Nn)2 - (^Ar)i]| as a function 
of the charge number and the neutron number of fission- 
fragments, we investigate the fissioning nuclei. For an 
example. Fig. Ufa) shows the neutron number versus the 
charge number of the fissioning nucleus for both reac- 
tion systems just before the fission takes place. The 
lines represent the second order polynomial fits to guide 
the eyes. From the above points, we can extract the 
[(Nz)2 — {Nz)i]\iis as a function of Zns as shown in 
Fig. El^b). Since the fission-fragments are assumed to 
have the same N/Z as one of the fissioning nuclei, the 
fission-fragments shall show the similar increasing behav- 
ior as the charge number rises. In comparison to the 
insensitivity of az or ctat to aao, [i^z)2 — iNz)i] and 
[(Nn)2 — {Nj^)i] shows somehow stronger dependence 



on aag, i.e. [{Nz)2 — {Nz)i] shows larger values in the 
proton-rich side while |[(Afjv)2 ~ (-^Af)i]| shows smaller 
values in the neutron-deficit side as becomes smaller 
in Fig. El This essentially originates from the overlap 
of the different mass partition of two fission-fragments 
according to the sampling of Gaussian probability dis- 
tribution for the mass asymmetry when the fission takes 
place. The symmetric fissions result in the strongest de- 
pendence of [{Nz)2 - {Nz)i] or \ [{Nn)2 - {Nn)i] \ on the 
charge number or neutron number. 



B. Isoscaling behavior 

Eq.(l) can be written as lni?2i = Cz + ctN, where 
Cz = InC -I- f3Z, if we plot R21 as a function of N, on a 
natural logarithmic plot, the ratio follows along a straight 
line. In Fig. [S] this isoscaling behavior is observed in the 
Langevin simulation. Here each kind of symbol with a 
line represents a chain of isotope. From there, the isoscal- 
ing parameter a can be extracted directly. Similarly, the 
isoscaling parameter (3 can be extracted from the isotonic 
ratio as shown in Fig. by lni?2i = Cn + (3Z, where 
Cn = InC + aN. 

From Fig. |S1 and Fig. |51 the relationship between a 
(|/3|) and the charge number Z (N) of the fission frag- 
ments can be deduced. In order to investigate the effect 
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FIG. 7: (Color online) (a) The calculated neutron numbers as 
a function of the charge numbers for the fissioning nuclei of two 
reaction systems, (b) The difference of the neutron number of the 
fissioning nuclei between two reaction systems as a function of the 
charge number for the fissioning nuclei. The lines represent two- 
order Polynomial fits in (a) and a linear fit in (b). 




FIG. 8: (Color online) The isotopic yield ratio of the fission frag- 
ments between ^^"^Sn -|- ^^^Sn and ^^-^Sn -|- ^^-^Sn in the Langevin 
model with = 0.06 and -Ebeam/A = 8.4 MeV. Different symbols 
from left to right represent the calculated results for the isotopes 
from Z = 37 to 59. The lines represent exponential fits to guide 
the eye. 
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FIG. 9: (Color online) The isotonic yield ratio of the fission frag- 
ments between -'^'^^Sn + -"^^^Sn and -'^'^^Sn -|- ^-"^^Sn in the Langevin 
simulations at Eheam 

I A = 8.4 MeV. Different symbols from left to 
right represent the calculated results for the isotopes from Af = 46 
- 73. The width of Gaussian probability is assumed 0.06. The lines 
represent exponential fits to guide the eye. 



of the width of Gaussian probabihty distribution on the 
isoscahng parameters, we use the different widths of the 
samphng Gaussian distribution for mass asymmetry pa- 
rameter ao, namely — 0.04, 0.06, 0.08 and 0.20, with 
the random number from -1 to 1 and the mean value of 0. 
Fig lior al shows the isoscaling parameter a as a function 
of Z with different ■ From this figure, we know in the 
low (Tq,q , i.e., when the symmetric fission is an overwhelm- 
ing mechanism, a increases with Z. This means that the 
isospin effect becomes stronger with the increasing of Z . 
In a recent analyse of fission with a simple liquid-drop 
model , a systematic increase of the isoscaling param- 
eter a with the proton number of the fragment element 
has been predicted. In our simulation, this kind of in- 
crease of a with Z apparently stems from the dominated 
symmetric fission mechanism. While, in the another ex- 
treme case from the Fig. Iior a). i.e. with the larger crc<oi 
a shows a contrary trend with Z, i.e., it drops with Z. 
In this case, it seems that there exists stronger isospin 
effect for the fragments with lower Z. In a medium 
case, the rising branch and falling branch competes with 
each other, the mediate isoscaling behavior appears and 
a minimum of a parameter occurs around the symmetric 
fission point. We note that the fission data of 238, 233^ 
targets induced by 14 MeV neutrons reveal the back- 
bending behavior of the isoscaling parameter a around 
the symmetric fission point as stated above. They in- 
terpreted that it originates from the temperature differ- 
ence of fission fragments since the isoscaling parameter 
is typically, within the grand-canonical approximation. 
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considered inversely proportional to the temperature ( a 
— Afin/T ) as stated above. In our case, this kind of 
backbending of isoscaling parameter a apparently stems 
from the moderate width of the probability distribution 
of the mass asymmetrical parameter of the fissioning nu- 
cleus as shown in Fig llOl or in the other words, it may 
stem from a moderate mixture of the different weights 
between the symmetric and asymmetric fission compo- 
nents. Essentially the backbending originates from the 
competition between the term of chemical potential and 
the term of temperature since both terms increase with 
the charge number of fission fragments. 

Besides the above direct method to extract isoscaling 
parameter, we can also check the behavior of a in terms 
of Eq. IISI). Fig.IUlb) shows the f(^^)2-(Afz)i1 ^ 
tion of Z. With the increasing of cTag, the Z dependence 
of shows from the upswing trend to down- 

swing trend. A turning point around Z — 51 is also 
observered in medium as Fig. llOf a) shows. From the 
similarity of the behavior shown in Fig. lior a) and (b) as 
well as the approximate equality of the values of a and 
[(Nz)2-{Nz)i] ^ gg^y ^^^^ ^j^g Eq.llll) works weU in 

the present calculation. In our case, the turning point of 
a stems from the competition between the chemical po- 
tential term {[{Nz)2 ~ (^z)i) and the temperature term 
(ct^). In general, the chemical potential term is more 
sensitive to the Gaussian width of the mass asymmetry 
parameter ao for fission fragments (see Fig. (B)). Overall 
speaking, we find that the isoscaling parameter a is sen- 
sitive to the width of the probability distribution of mass 
asymmetrical parameter of the fission fragments. In the 
other word, we may say that the isoscaling parameter 
is sensitive to asymmetrical extent of both fission frag- 
ments. 

Similarly, from Figure. |H| the relationship between |/?| 
and the neutron number N of the fission fragments can 
be deduced in different width . This is shown in 
Fig. Ilir a). Different from the relationship of a and Z, 
the always drops with the neutron number, regardless 
of the change of ■ The quantitative and qualitative 
similarity of ^ (Figure. [lirb)) has also 

been observed, i.e., it always decreases with TV and is 
insensitive to aa^. 

However, the obtained isoscaling parameters are actu- 
ally very large in comparison to the usual isoscaling pa- 
rameter extracted from the data 0- The reasons could 
be the model itself since the model is still too simple as 
well as our assumption of Gaussian probability distribu- 
tion of fission fragments. Also the post-fission evapora- 
tion component will of course play some roles for modi- 
fication the isoscaling parameters. In the present model 
calculation, however, this influence of post-fission evap- 
oration of fission fragment is not included. Those may 
show larger apparent isoscaling parameters in compari- 
son to the data. Of course, main aim of this work is to 
show the isoscaling behavior of fission fragments and its 
trend with the charge or neutron number of the fragments 
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FIG. 10: (Color online) (a) The isoscaling parameter a as a func- 
tion of Z in the different Gaussian width ((Tcq ) of the mass asym- 
metry parameter ao for fission fragments; (b) Same as (a) but for 

[(JVz)2-(JVz)ll 



by the Langevin dynamics. 



C. The beam energy dependence of the isoscaling 
parameters 

The simulations are systematically performed in differ- 
ent beam energies. The values of a and (3 are extracted 
as a function of beam energy for the fragments Z — - 
54 and = 58 - 68, respectively, as shown in Fig^Ja) 
and (b). It shows that both a and (3 decrease as the in- 
creasing beam energy which means that the isospin effect 
fades away with the increasing of Ehcam/A. This behav- 
ior is similar to the case in the fragmentation where the 
isoscaling parameter drops with the temperature in the 
statistical models as well as experiments |^ 0, |3^ 0| . 



D. The friction parameter dependence of the 
isoscaling parameters 

In addition, the influence of the reduced friction pa- 
rameter on the isoscaling parameters is investigated, we 
use a constant value of f3o = 2, 4, 6, 8 and 10 instead of 
one-body dissipation Pqbb which was used in above cal- 
culations. In Fig. Il^f a) and (b), we plot a and \(3\ as a 
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FIG. 11: (Color online) (a) The isoscaling parameter as a 

function of Z in the different Gaussian width (itcio ) of mass 

asymmetry parameter ao for fission fragments; (b) Same as (a) 
but for (^")2-(iViv)i |^ 



FIG. 12: (Color online) Isoscaling parameter a (a) and fi (b) as a 
function of beam energy for the fragments Z = 44 - 54 and A'^ = 
58 - 68, respectively. The width (0^^ ) of the Gaussian probability 
is 0.06. 



function of /3o for different elements from Z = 44 to 54 or 
different isotones from N = 58 to 68, respectively. Both 
a and |/3| decrease with the increasing of the reduced fric- 
tion parameter. It shows that a and /3 are sensitive to the 
the reduced friction parameter. Larger reduced friction 
makes the Brownian particles cost more energies which 
will be transferred to the internal energy from ground 
state to the scission point than the smaller one, conse- 
quently the system will keep less memory to the initial 
entrance channel. In the viewpoint of the isoscaling be- 
havior, the isoscaling parameter shows a decrease with 
the reduced friction parameter. Therefore the study on 
the isoscaling behavior to the fission fragment might be 
a good tool to explore the friction effect in the fission 
dynamics process. 



IV. SUMMARY 

In summary, we applied the Langevin model to inves- 
tigate the isoscaling behavior in the dynamical process 
of compound nuclear fission. In order to treat the fission 
fragments, we assume that the mass asymmetry parame- 
ter of the two fission fragments from the fissioning nucleus 
is taken from a random number with a Gaussian distribu- 
tion whose width is Uan- The simulation illustrates that 



the isotopic and isotonic yield ratios of fission fragments 
in the dynamical fission channels of ^^^Sn -|- ^^^Sn and 
^^^Sn -I- ^^^Sn reaction system show the isoscaling behav- 
ior. The terms which are related to the difference of neu- 
tron or proton chemical potential are also extracted. It is 
of interesting that the isoscaling parameter a is sensitive 
strongly to the Gaussian width CTq^ of the mass asym- 
metry parameter but /3 looks not. When CTa^ is small, 
i.e. the fission is almost symmetric, a increases with the 
atomic number of fission fragments, which is similar to 
the theoretical prediction of a simple liquid-drop model 
jsj. In contrary, when is large, for instance, = 
0.20, a drops with Z of fission fragments. However, in 
the intermediate values of ctq,-,, a shows a backbcnding 
with Z of fission fragments, which is similar to the ob- 
servation of the 238,233y figgion data induced by 14 MeV 
neutrons 0. In this context, we could say that the a 
parameter is sensitive to the asymmetric extent of the 
fission-fragments from the fissioning nuclei. However, fj 
parameter is insensitive to the width even though it 
always shows the dropping trend with N . 

In addition, the dependences of beam energy and the 
reduced friction parameter for the isoscaling parame- 
ters are systematically investigated. It is found that 
both a and /3 drop with beam energy of the projectile 
as well as the reduced friction parameter, refiecting the 
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FIG. 13: (Color online) a and \0\ as a function of the reduced 
friction parameter I3q for different fission fragments with Z = AA - 
54 and = 58 - 68, respectively. The width {o-ag ) of the Gaussian 
probability is 0.06. 



temperature-like dependence of isoscaling parameters in 
the fission dynamics. The disappearance of isospin effect 
of fission dynamics is expected in a certain higher beam 
energy or larger reduced friction parameter. In general, 
the isoscaling analysis of the fission data appears to be a 
sensitive tool to investigate the fission dynamics. 
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